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ABSTRACT 



An impcoved and simplified modal £oc the 
Heave-Pitch dynamics of the XR-3 CA3 Surface Effect 
Ship is developed for constant speed operation. The 
nonlinear equations of notion ace linearized aoout the 
steady-state operating point. Time-domain validation 
is accomplished by comparison with the 6 DO? model. A 
signal flow graph of the craft dynamics is developed 
and ilason's Sain Rule used to determine the 
characteristic S-Domain polynomials describing the 
craft's vertical plane dynamic behavior, with Bode 
Plot analysis included. Conclusions are drawn and 
recommendations for further study are made. 
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I. 



INTRODUCTION 



In 1973 Leo and Boncal [1] converted the Oceanics Inc. 6 
DOF Loads and Motions digital computer simulation program 
for the 100-B Surface Effect Ship to reflect the physical 
characteristics of the XR-3 craft. ' Since that time there 
has been a continuous effort at the NPS to update and 
improve on this simulation [2, 3,4,5] . 

One of the major drawbacks of that program which 
precludes regular use of it for sea stats operations 
simulation is the excessive computational and turn-around 
times required [6] . As it is presently desired to use the 6 
DOF simulation for design and evaluation studies of the 
proposed 3K Ton 5SS control System in varying sea state 
conditions, it is important to significantly reduce the 
simulation's computational time. 

One area where the possibility exists to achieve a major 
reduction in computational time is the manner in which 
sidewall forces and moments are calculated. The present 
form of the 6 DOF program uses a large number of discrete 
sections (28 per sidewall) with an off-line table look-up 
scheme to determine the forces generated in each discrete 
section. The simple model derived by McIntyre DJ uses two 
sections per sidewall with very simplified geometry that 
permits on-line computation of the forces and moments. The 
result of eventual incorporation of this sidewalL geometry 
and forces calculation into the 5 DOF program should be a 
significant reduction of computation and turn-around time. 

Another area of great interest in SES Simulation is 
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gaining insight into the dynamic behavior of the craft's 
vertical motion of pitch and heave due to pressure induced 
lift forces generated in the plenum and seals. is stated in 
ref C8] as recently as August, 1 978, 

'•The overall precision of the analysis is satisfactory 
for predicting stability and maneuverability. However, 
several areas present very complex or difficult physical 
phenomena to analyze and further effort should be 
invested to improve understanding of these. ' ' 

The first area for continuing investigation recommended was; 

•'A better determination of pitch-heave dynamic 
characteristics. ' ' 

An improved knowledge of the dynamic characteristics 
brougat about by the pressure induced forces and moments 
could lead to a reduction of the computation time 
required for craft simulation. These pressure effects 
introduce a stiffening in the integration of sea state 
induced oscillations and a better understanding of 
processes could lead to the removal of this stiffness 
and thus reduced CP r J time. 

In this work both of these areas are considered. First 
of all the various craft physical characteristics are 
modeled, resulting in changes and improvements to the 
initial efforts by McIntyre C7J . These equations are 
linearized about a steady-state operating point using 
the Taylor Series expansion. The linearized equations 
are cast into a State Space format and simulated on the 
digital computer. Comparison is made between the 

time-domain results of the linear model and the 6 DOF 
model. The State Space representation is used to 

generate a signal flow graph. Mason's Gain Rule is 
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applied and the S-domain roots of the craft's 
characteristic equations are located and it is shown 
that some simplification and reduction of system order 
is possible. The resultant S-domain polynomials are 
used to generate Bode Plots to graphically display the 
craft's simplified frequency response. 

Conclusions are drawn and recommendations are made 
concerning the validty of the model and areas for future 
study . 

Several procedural notes are in order. Because all 
equations are to be solved digitally using the FORTRAN 
computer language and because digital computer text 
processing facilities were used to generate this 
document, all formulas will be presented in FORTRAN 
format. In the thesis body reference to the ''nonlinear 
modal'', ''six degree-of -freedom model'' or ''6 DOF'' 
model all mean the Loads and Motions program as adapted 
to the XR-3 test craft at the Naval Postgraduate School 
C3J . Reference to the linear model refers to the 
linearized modal of the same craft as developed herein. 
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II. NONLINEAR XR-3 MODEL HEAVE AND PITC H EQUATIONS 



A. BACKGROUND 



The original development work on the simplified 
two-degrae-of -freedom model was done by McIntyre [7] .. 
Several weaknesses in the original model were noted by 
McIntyre in that work. 

The major portion of the work presented here was 
involved in close examination of all aspects of the 
original model, comparison of the modeling techniques 
used therin with other methods and remodeling of several 
important craft physicaL characteristics. 

One major area of investigation was the lack of pitch 
damping in the original model, which resulted in the 
addition of an added-mass effect to ' the pitca moment 
equations . 

Additional changes involve a revision of the force 
aquations for the ceac seal, major changes la the way 
the plenum chamoer is modeled (also added to the loads 
and motion 6-dof model C23 | . and a more rigorous 
development of the planing force equation of [7] . 

There are several areas of theory concerning CAB 5ES 
modeling where empirical data is lacking due to 
difficulties with physical measurements. One such area 
inviting closer inspection was dynamic variation of the 
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effective canter of pressure for the plenum gauge 
pressure lift effect. 1 possible approach is suggested 
herein and included in the model. 

B. ASSUa?TIONS 

The following assumptions are made for this constant 
speed, pitch-heave modal: 

1. Planing forces ara concantratad at tha cantroid of 
the keel longitudinal cross section. 

2. Total saal force is lumped into a pressure 
differential acting on the seal face area in contact 
with the water passing baneath the saal. 

3. Aerodynamic forces on tha forward seal are 
disregarded. 

4. Dynamic variations in momant lever arms are 
disregarded for all monants except seals and tha moment 
due to plenum guage pressure. 

5. Dynamic variations in effective plenum roof area 
will be included. 



6. Any lift effects 


due to 


aerodynamic 


forces are 


lumped in with the 


planing 


forces 


{for 


reasons to be 


developed in Chapt IV. 


) The 


moments 


a re 


included in 


plenum gauge pressure 


moments 


• 







7. Constant seal leakage area is assumed. 

8. Since this model is for constant surge velocity. 
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drag and thrust forces are considered to be in 
equilibrium and in balance for moments about the Y axis. 



C. COORDINATE SYSTEM 

(see fig. 1 ) 

1. The origin of the cartesian coordinate system is 
located at the calm-water line, longitudinally and 
transversly coincident with the craft center of gravity. 

2. The X axis is the longitudinal coordinate defined 
as positive as X increases towards the bow of the craft. 

2. The Y axis is the transverse coordinate increasing 
positively towards the starboard side. 

4. The Z axis is the vertical axis, increasing 
positively in a downward direction. 

5. ZS is the vertical distance from the center of 
gravity to the keel. 

5. Z is the vertical distance from the calm water line 
to the vertical canter of gravity. 

7. THETA is the pitch angle defined as positive as the 
bow pitches up. 

8. Positive pitch moments are defined as moments which 
tend co cause a positive increase in the angle THETA. 



D. SIDEWALL MODELING 
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As it was desired to keep ail physical dimensions 
consistant with the six degree of freedom model, scale 
model drawings were developed with the data from the 
sidewall subroutine of that modal. All sidewall 
dimensions used in this model were taken from those 
drawings. 

The sidewalls are modeled in two sections, one forward 
of the longitudinal 3.3. and the other aft. Each 
sidewall is considered to be of constant width at the 
water line and the keel, with width increasing in the 
vertical direction according to an average deadrise 
angle. 

The average forward sidewall area is found to be: 

(see fig. 2 and 3b ) 

AVERAGE WIDTH = LDB AR1 / (2 * TAN (DR 1 ) ) <• WS 1 0 



(II-1) 

where LDBAR1 is draft at forward centroid, DR 1 = 

average deadrise angle forward and WS10 is the average 
keel width of the forward section. Therefore submerged 
volume of a forward sidewall is given by: 

VOL = LI* (LDBAR 1/(2*TAN (DR1) ) +WS10) *LDBAR1 



(H-2) 



where LI is the length of a forward buoyant section. 
Following the same procedures the after sidewall width 
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is: 

(see fig. 2 and 3c ) 

AVERAGE WIDTH = LDBAR2/ ( 2 *TAN (DR2)) * WS 2 0 

(II-3) 



where DR2 = average deadrise angle aft 
LD3AR2 = draft at after buoyant centroid 
WS20 = average Iceel width aft 

Submerged volume of an after sidewall is then 
determined to be; 

VOL = L2* (WS20 + LD3AR2/(2*TAN(DR2) ) > *LD3AR2 



(H-4) 



where L2 is the length of the after sidewall. 



E. SEAL MODELING 



A major change in McIntyre's model was made by 
remodeling the stern seal. In the original model a large 
negative moment was required to place the pitch 
equations in steady state equibibrium. Investigation of 
seal forces as presented in C9Q revealed a large 
discrepancy between the after seal forces generated in 
McIntyre's model and calculations based on the 

information in BO . This discrepancy was further 
supported by actual seal load data presented by Layton 
in Cl CD . All seal dimensions used in this model were 
taken from data presented in D "0 . Seal hinge locations 

were taken from the data used in G] . 



17 



Seal lift force forward is generated by plenum gauge 
pressure acting downward on the seal face area in 
contact with the water passing underneath the seal. 
This downward force results in a lift force being 
generated as the downward movement of the seal displaces 
the water beneath it. 

The seal is considered to have a straight forward face 
which runs from the hinge point -at the plenum roof to 
the Qottom of the keel at an angle determined to be 
best-fit from the data in ref. DTJ . The wetted surface 
of the seal is considered to lay parallel to and on the 
water at the point where the seal face intercepts the 
calm water line, as a result forward seal lift equations 
are : 



ASEAL1 = WIDTH*XSEAL1 



(II-5) 



XSSAL1 = (LD-L3*TAM (THETA) ) /SIN (31) 



(II-6) 

LD = draft at the C.G 

L3 = distance from Z.G to point of contact for forward 
seal face. 

LD - L3*TAN (TKETA) = draft at forward seal face contact 
point 

31 degrees is the interior angle of the seal. 

The after seal watted area is determined in a similar 
manner to be : 

AS SAL 2 = WIDTH* XSEAL 2 
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(H-7) 



XSEAL2 = (LD+L4*TAN (THETA) ) /SIN (32) 



(I I- 8) 



In this model forward seal pressure is considered 
equal to plenum pressure. As a result forward seal lift 
is : 



HSEALF = PB3 AH* ASEAL 1 



(II-9) 



Since plenum guage pressure tends to force the after 
seal upward, resulting in venting of the plenum air 
mass, the rear seal is operated at a pressure slightly 
above that of the plenum to reduce the venting tendency 
of the after seal. As a result the after seal requires a 
modification in modeling. In McIntyre's work the 
pressure differential (PSE AL-PBBAR) only was considered 
in the lift force. However, as previously stated, there 
were large discrepancies in computed loads as compared 
to measured data from Cl OH . Therefore in this model , 
the seal working pressure is considered to be P33AR plus 
differential pressure. As a result after seal lift is 
given by: 



HSEALA = (PBBAR + PDIFF) *A5EAL2 



( 11 - 10 ) 
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F. PLENUM MODELING 

In the work presented by McIntyre the plenum was 
considered a box with the only dynamic dimensional 
variation being in the Z direction. However (as can be 
seen in fig. 1 ) , although the forward edge can be 

assumed to be reasonably constant in the Y-Z plane, the 
forward face of the rear seal slopes rather steeply. 
Therefore the plenum length (in the X direction) varies 
rapidly as draft varies. As a result the effective 
plenum roof area consists of that area from the after 
edge of the forward seal to the forward hinge point of 
the rear seal, plus the vertical projection of the 
unwetted forward face of the rear seal. As an 

additional consequence it can be seen that the effective 
plenum center of pressure varies fore and aft as draft 
varies . In fact as draft increases the Z.P. moves 
forward. 

The following variable names will be used to describe 
the plenum chamber (see fig. 4.) 

LP8 = length of plenum roof from the rear eije of the 
forward seal to the hinge point of the after seal. 

LPWL = length of plenum at rhe keel from the after edge 
of the forward seals to the transom. 

LP = length of plenum at the waterline. 

3U3HGT = vertical distance from the plenum roof to the 
bottom of the keel at the center of gravity. 

A3 = plenum effective roof area. 

Using the previous definitions, the plenum roof length 
and area is: 
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LP = LPWL- (LPWL-LPR) *(LD/BUBHGT) 



(II- 11a) 



and effective plenum roof area is 
AB = WIDTH *LP 



(II- 11 b) 

with the result that lift due to planum gauge pressure 
is: 

HPRES = -A3*PBBAR 

(II-llc) 

As a comment on the complexity of the lift forces of 
this model note that for constant plenum pressure 
increased draft results in decreased plenum lift. 

G. CRAFT DYNAMICS 



1 • A ir Mass and Pressure Fguat ions 

For the purposes of this model the rear seal 
leakage area is considered constant since to assume 
otherwise would lead to modeling a very difficult and 
poorly understood dynamic process. In addition the 
plenum pressure VS air mass process is assumed to be 
adiabatic in nature. Fan characteristics are tiose used 
in the constant air leakage heave only model of Gerba 
and Thaler in C143 . Air flow into the plenum is taken 
as positive. 
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Therefore the airflow into the chamber is given 

by: 

Qin = N*(QlO-Kq*PBBAR) 

( 11 - 12 ) 



Where N is the number of fans supplying air to the 
plenum, QIO is the equilibrium air flow rate for the fan 
with zero gauge pressure and Xq is the slope of the 
fans' airflow VS gauge pressure curve . 

Qout is airflow out of the plenum and is 
determined oy the seal Leakage equation: 

Qout = Cn*Al* (2*P5BAR/RH0A) **1/2 



(H-13) 

Where Cn = seal leakage coefficient, Al = seal leakage 
area, and RHOA = mass density of air. Therefore mass 
flow rate is: 

MB DOT = RHOA* (Qin-Qout) 



(11-14) 



The plenum air dynamics are considered in 
equilibrium when Qin = 2out 

Absolute Plenum Pressure is determined from the 
adiabatic gas law to be: 

Pb = Pa* (MB/ (RH0A*PB3AR) ) **GAMHA 
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(H-15) 



and. 



PB3AH = Pb-Pa 



(II- 15a) 

Where Pa = ambient atmospheric pressure (LBF/Fr**2 GAMMA 
= adiabatic process coefficient, Pb = plenum absolute 
pressure, and PB BAR = plenum gauge pressure. 

2 • Buoyant Fo rces 

Buoyant force is calculated by the submerged 
volume multiplied by the mass density of the fluid times 
the gravitational constant. 

Therefore forward buoyant force is: 

HBF = -2*RH0*G*V0L 



(11-16) 

Where volume is given by eqn. II-2, RIO is the 
mass density of water (slugs/f t**3 ) and 3 is the 

gravitational constant. The factor of 2 is introduced 
to account for 2 sidewalls. The total expression for 
HBF is: 

HBF = -2*RH0*G*L1*LDBAR1 * (LDB AR1 / (2*TAN (DR1) ) +WS01) 

(11-17) 



Following the same pattern. 
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HBA = -2*RH0*G*L2*LDBAR2* (LDBAR2/(2*TAN (DR2) ) =W51 0) 

(11-18) 

3 . Seal Forces 

Forward seal force is determined from eqn. II-9 to 

be : 

US F = -P3BAR*WIDTH* (LD-L3*T'AN (THETA) ) /SIN (31) 

(11-19) 

where HSF is the heave force due to forward seal. 

In a similar manner HSA is: 

HS A = (PB8 AR + PDIFF) * WIDTH* (LD+L4*TAN (THETA) /SIN (32) 

( 11 - 20 ) 



4 . Pla n ing Forces 



McIntyre introduced the planing force phenomena in 
his model to account for the proven fact that the CAB 
SES tends to decrease it's draft and pitch down as speed 
increases. His development of the planing equation was 
not well explained and the basis for selection of the 
planing coefficient was not clearly stated. Reidel Cl 2] 
did a special investigation of the planing force 
phenomena. His results were based on flat plate planing 
effects to approximate the flat bottom of the craft 
sidewalls. 
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Since the water is in contact with the vertical 
portions of the sidewalls it was felt a more accurate 
representation was neeced. Barnaby Cl 33 states that 
water contact on the vertical sides of a planing hull 
rapidly reduces planing lift advantage since one of the 
advantages of a planing hull is reduced drag due to 
wetted surface. planing lift results in a speed 
sensitive tradeoff between static lift (buoyant effects) 
and dynamic lift (planing effect) that is directly 
proportional to the angle of attack (for small planing 
angles) and the square of the velocity. Therefore the 
planing force is of the form: 

HPLAN = K*V**2*THETA 

From Reidel's work, 

HPLAN = 0. 5*AREA*V**2*PI*SIN(?HETA) 



( 11 - 21 ) 



for small angle approximations this simplifies to: 

HPLAN = RHO*AEEA*V**2*?I*THETA 

Which is the planing force for a flat plate. Barnaby 
gives no simple formula for the losses due to immersion 
so we are left with an arbitrary loss coefficient, 
PLCOEF such that: 

HPLAN = (PLCOEF) *2 *RHO* AR 5 A* V**2*P I*THET A 

( 11 - 22 ) 

The added factor of 2 accounts for 2 sidewalls. 

The planing area is taken from the scale drawings 
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developed from sidewall data, and the force is assumed 
to act at the centroid of the flat keel surface, which 
is consistent with Reidel's work. The loss coefficient 
will be developed in the section on achieving 
equilibrium. 

5. Ple n um Gau<qe Pressure Lift 



Combining equations Il-lla ,II-11b , and II-llc 

results in: 

HPRES = -PBBAR*WIDTH*(L?WL- (LPWL-LPR) * (LD/BUBHGT) ) 



(11-23) 



6 . Total Lift For ces 

Lift Total = HPRES + H3F +• HBA + HSF + HS A ♦ HPLAN 

which must exactly equal craft weight (W) for the hull 
to be in vertical equilibrium, or; 

W *■ HPRES + HBF + HBA + HSF + HS A + HPLAN = 0 

(11-24) 

defines equilibrium condition in heave motion. 



H. CRAFT PITCH MOMENT EQUATIONS 



1 • Buo y ant Mome nts 
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Forward buoyant moment is the forward buoyant lift 
times the effective lever arm, or; 



PBF = -HBF* (Lever arm forward) 



PBF = -HBF*L5 



(11-25) 



Note that the 
definition of positive 
forces are negative by 



negative sign results from the 
moments and the fact that upward 
definition. 



PBA = HBA*LEVER AR ii 



P3A = HBA *L6 



using the same reasoning. 



(11-26) 

after buoyant moment is: 



PB total = PBF + PBA 



(11-27) 



2 • Se al Moments 



C.G 

in 



Since the seals ace relatively distant from the 
. dynamic variations in the lever arms are considered 
seal moment equations and are developed uere. 
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The seal force is modeled to act at the 
longitudinal center ofthe seal wetted length. Therefore 
the effective lever arm for the seal moment is: 

PLSF = L3 - XSEAL1/2 

where PLSF = moment lever arm for seal forward. As a 
result : 



PSF = -HSF*PLSF 



PSF = -HSF* (L3- (XSEAL1/2) ) 



(11-28) 



In a similar manner 



PLSA = L4 + XSEAL2/2 



PSA = HSA*?LSA 



PSA = HSA* (L4 + XSEAL2/2) 



(11-29) 



PS TOTAL = PSF+PSA 

Note that the dynamic pitch sensitivity to moment 
arm length is included in the XSEAL terms as defined in 
equations II-S and II-9. 
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3 • Planing Moment 

For the planing force developed in equation 11-22 
the planing force is concentrated at the centroid of the 
flat keel which results in: 

PPLAN = L7*HPLAN 



(11-30) 



4. Plenum Pressure Moment 

First the dynamic variation in affective center of 
pressure will be developed, (see fig. (4)) 

plenum length is taken from aguation Il-lla to be 

LP = LPWL- (LPWL-LPH) * (LD/BUBHGT) 

Defining LDIFF to be 

LDIFF = LPWL-LPS 



(11-31) 



and. 



XCP = XCP0+ (LDIFF*LD/3UBHGT) /2 



where 



(H-32) 
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XCPO = L PWL/2-XCG 



(11-33) 

and XCG is the distance from the transom to the C.G. 
(see fig. 4 .) Finally, 



PPRES = - HPRES*XCP 



(11-34) 



5. Pitch Damping Mo ment 

Any frictional damping terms dependent upon the 
square of pitch rate (i. e. tangential velocity about the 
I axis) disappear when the Taylor Series expansion is 
carried about the zero rate equilibrium point. Since 
for small signal approximation the model should reflect 
actual craft characteristics, it was hoped that the 
changes in seal and plenum modeling would reveal the 
missing damping characteristics. As it turnei out this 
was not the case. Examination of the hydrodynamic 
equations in GT] revealed an added mass effect dependant 
upon a cross velocity term between linear X direction 
velocity U, pitch rate Q, and added-mass at the stern 
A33S due to slender body theory applied to the sidewalls 
. This moment is: 

FP = -A33S*XLSS**2*U*T3ETADOT 

(11-35) 

Where XLSS is the distance from the C.G. to the stern, 
and TH2TAD0T is pitch cate. 
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A3 3S = YSS**2*RHO*PI/8 



(11-36) 

where YSS is the :ross section width of the side 
wall at the stern waterline. Therefore, 

FPDAMP = -2*XLSS**2*A33S*U*THETADOT 



= - 2*PI*XL5S**2*YLSS**2*RH0*U*THEID3T /8 

(11-37) 

The factor of 2 is introduced to account for two 
sidewalls . 



6 . Total Mome nts 



Since the moments are all referenced to the craft 
center of gravity, the craft weight contributes zero 
moment in that frame of reference. Therefore total 
moments are: 

PSF + PSA + P3F + P3 A * PRES + PPLAN + PDAMP = 3 

(11-38) 



I. HEAVE ACCELERATION EQUATIONS. 



From Newton's law, acceleration and force are related 
by: 
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FORCE 



MASS* ACCELERATION 



For the linearized model we are interests! in the 
acceleration on the craft and as a resalt: 



ACCELERATION = FORCE/ MASS and therefore the heave 
equation of motion becomes: 



ZDDOT = (HBF+HBA+HSA+HPL AN+HPRES) /MASS 

where ZDDOT is the second derivative of 
to time. 



( 11 - 39 ) 
Z vita respect 



J. PITCH ACCELERATION EQUATIONS 

Angular acceleration is: 

THETADDOT = PTOTAL/IYY 

Where IYY is the craft inertial moment about tae Y axis. 
As a result total angular acceleration is; 

THETADDOT = (P3F *PBA +PSF + PS A+ PPL AN + P PR ES +P DAMP) / IYY 

( 11 - 40 ) 



This completes the derivation of all nonlinear 
equations needed to describe the two-degree-of- freedom 
CAB SES. In the following chapter these nonlinear 
equations will be linearized using the Taylor Series 
expansion about a steady state operating point which 
will result in the equations needed to generate a 
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notion in 



dynamic digital-computer simulation of craft 
the time-domain. 
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III. LINEARIZATION OF THE HEAVE-PI TCH 32UA££ONS. 



A. TAYLOR SERIES EXPANSION 



For small perturbations about a steady state operating 
point the nonlinear function Z can be approximated by: 

Z = Z (X,Y) 

Z (0) +dZ = Z (X (0) ,Y(0) ) f ( JZ/iX) *dZ* (dZ/dY) *dl 

where ($Z/dX) and (&Z/dY) are the partial derivatives 
of the function Z evaluated at the steady state 
operating point X(0),Y(0). Cancelling the steady state 
values from both sides of the equation yields the 
differential model equations: 

dZ = (&Z/dX) *dX+ (AZ/ZY) *dY 
and 



Z = Z (0) +dZ 

where dZ r dX r and dY are the differential variables for 
Z, X, Y respectively, (3Z/<}X) and (AZ/3Y) are the 
partial derivatives of Z evaluated at the steady-state 
operating point. 
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B. LINEARIZATION OF TBS MODEL EQUATIONS 



There are basically four equations to linearize for 
the two-degree-of- freedom model: 

ZDDOT = F (Z,THEIA, MB) 



THSTADDOT = G (Z , THETA, T3 ET ADOT , MB) 



MBDOT = 3 (Z r MB) 



P3BAE = I (Z , MB) 

therefore the equations for the differential state 
variables are: 

d ZD DOT = (dF/dZ) *dZ + (dF/dTHETA) *dTHETA+ (dF/dM3) *dM3 

(III-1) 

dTHETADOT = (dG/dZ) *dZ«- (dG/dTHETA) *dTHETA+ (dG/dM3) *dMB 

(HI - 2) 

dM3DOT = (dH/dZ) *dZ + (dU/dTHETA) *dT r dETA + (dH/dM3| *dMB 

(III-3) 
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dP3BAR = (dl/dZ) *dZ+ (dI/dTH3TA) *dIHETA+ (dl/daB) f'dMB 

(III-4) 

Mote that the left hand side variables are 
differentials, not time derivatives. 

C. LINEARIZED EQUATIONS 



1. Plan aim A ir Mass F L ow Derivativ es 
Qin = N*(QiO-Kq*P33AR) 

( 11 - 12 ) 

Qout = Cn*Al* (2*PBBAR/RH0A) **1/2 

(H-13) 

MBDOT = a HO A* (Qin-Qout) 

(11-14) 

as a result the differential equation for MBDOT becomes: 
dMBDOT = RHC A* (dQin-dQout) 

(HI-5) 

iQin = (dQin/dZ) *dZ + (dQin/dPBBAR) *dP3BAR+ndQin/dMB) *dMB 
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(Ill - 6) 



dQin/dPBBAR = -N*Kq 

The partial derivatives of Qin with respect to all other 
variables are zero. 



dQout = 

(dQout/dZ) *dZ+dQout/dP3BAR) *dPBBAR+ (dQout/dMB) *dMB 



dQout /dPBBAR = Cn*Al*((3flOA/2*PB3AH) **1/2) /RHO A 



(III - !) 

The partial derivatives of Qout with respect to all 
other variables are zero. Inserting equations EII-6,and 
III-7 into 11-14 yields the equation for the 
differential variable ddBDOT: 

d M 3 DO T = (-RHOA*N*Kq-Ca*Al* ( (RHOA/2* PBB AR) **1/2) ) *dPBBAR 

' (H-15) 



2. Ple nun Gauge Pressure Derivatives . 



Pb = ?a*(db/RHOA*Vb) **3A?5tlA 



(III-8) 



and 



PBBAR = Pb-Pa 
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Since Pa is considered constant within the time frame of 



this model. 



dPBBAR = dPb 



(III-9) 



As a result. 



dPBBAR = 

GAdilA*Pa* (MB (0) /RHQA*Vb (0) ) ** (GAMMA- 1) *d (HB/RHOA* VB) 
where 

i (33/RHOA * V3) = (1 /RHOA* VB (0) ) *dMB+ (MB/RHOA) *d ( 1 /VB) 

d(33/RH0A*Vb) = (1/RHOA*73(0) ) *dM3- ( MB/RHO A*V3** 2) *dVB 

dV3 = d(VN-AS*LD) 



dV3 = -AB (0) *dLD/dZ-LD ( 0 ) *dA3/dZ 
remembering that: 



LD = Z + ZS 
then 



aLD/dZ = 1 



we now need the partial of affective plenum roof 
area with respect to Z. Defining 
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AB = ABrfL- (ABWL-ABR) *LD/BUBHGT 



(III — 1 0) 



ABWL = WIDTH*LPW L 



ABH = WIDTH*LPR 
farther defining 

ABDIFF = ABWL-ABR 

which is a constant, the effective plenum roof area is: 



A3 = ABWL-ABDIFF*LD/B(JBH3r 

using the previous definitions the partial of effective 
roof area becomes: 



dAB/dZ = -ABDIFF/3U3HGI 



and the differential volume is: 



dVb =- (AB (0) -LD (0) * 2 * ( A3 D I FF/BCJBHGT) ) *dZ 



with the result that 



(III— 1 1) 



(III-1 2) 



d (M3/RH0A* Vb) = (1/RHOA*Vb(0) ) *dZ * 



(Mb (0) /3HOA*Vb (0)**2)*(AB(0) -2*LD (0) * (ABDIFF/B J3HGT) ) *dZ 
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finally , 



dPBBAR = GAMMA*Pb(0) *((1/Mb(0) ) *dMb * 

(1/Vb(0) * (AB (0) -2*LD (0) * (ABDIFF/3 UBHGT) ) ) *dZ 

(III-13) 



3- Heav e Force Derivatives 



The total heave force equation is: 



H (total) = HBF+HBA + HSF+HSA+HPLAN+HPEES 



The total partial derivatives will be done term by 

term. 



For the forward buoyant force the following 
simplifications are made: 

H3F = -2*R HO*G*L 1* LD BARI * (LDBAR1/(2* TAN (DR 1) ) + W510) 
let K 1 = 2*RHO*G*L1 then, 

HBF = -K1*LDBAR1* (LDBAR1 / (2*TAN (DR1> ) +MS10) 

where LDBAR1 = (LD-L5TAN (THETA) ) . Combining the 
previous two equations, 

HBF = - K1 * ( LD-15 *TA N (THETA) ) * 

( (L D-L 5 *TAN (THETA) / (2 * T A N (DR1) ) f«SlO) 
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angles. 



assuming small pitcn 
approximated by THETA which 
simplification: 



TAN (THETA) can be 
results in the following 



HBF = 

-K1* (LD-L5*TH2TA) * ( ( LD-L5 *THET A) /(2*TAN (DB 1) ) +W3 10) 

Taking the partials of HBF with respect to Z and 
THETA yields: 

d (HBF) /dZ = -K1* ((LD (0) -L5*THETA(0) /TAN (DR1) + WS10 

(ITI-14) 



and 



d (HBF) /dTHETA = 

K1*L5* ((LD (0) -L5*THETA (0) ) / (2*TAN (D/rl) ) +HS10) 
+ K1* (-L5/(2*TAN (D31) ) * (LD(0) -L5*THET A (0 ) ) 

Which simplifies to : 



d (HBF) /dTHETA = 

K1*L5* ( (LD (0) -L5*TH2TA (0) ) / (TAN (DR 1 ) ) +WS 10) 



(III-15) 

Using the same approach, the partials for buoyancy 
force aft are: 

d (HBA) /dZ = -K 2* ( (LD (0) + L6*THETA (0) ) / ('TAN (DR2) ) t- WS20 ) 

(III- 1 6) 
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d (HBF) /dTHETA = 

-L6*K2* ({LD (0) + 16*1 H ETA (0) / (TAM (DR2) ) +WS10) 



( I II- 1 7) 

For the forward seal the total differentials are: 

DHSF = (d(HSF) /dZ) *dZ+ (d (HSF) /dTHETA) *d TH E T 
♦ (d<HSF) /dPBBAS) *dPBB AR 

the following simplification is made; 



HSF = ~P3BAR*NIDTH* (LD-L3*TAN (THETA) /SIN (3 1) 

(11-19) 

let WIDTH/SIN (31) = K3 and TAN (THETA) = THETA, then, 

HSF = -K3*?BBAR* (LD-L3*IHETA) 
with the resultant partials being: 



d (HSF) /dZ = -K3*PBBAR (0) 



(III-18) 



i (HSF) /dTHETA = K3*L3*PBB AR (0) 



(III-19) 



d (HSF) /dPBBAR = -K3* (LD (0) -L3*THETA (0) ) 



(III-20) 
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The partials of after seal force are found in the 
same manner to be: 

dHS A = (d (HSA) /dZ) *dZ+(d(HSA) /dTHETA) * dTHETA 
+ (d (HS A/dPBBAR) * dP B3 A R 



HSA = - (PBBAR+PDIFF) *WIDTH* (LD+L4*TAN (THETA) /SIN (32) 

( 11 - 20 ) 

let K4 = WIDTH/SIN (3 2) and assume small THETA so that; 



HSA = -K4* (PBBAR + PDIFF) * (LD+L4*THETA) 
the partials then become: 

d (HSA) /dZ = -K4* (PBBAR (0) +PDIFF) 

(III-21) 



d (HSA) /dTHSTA = -K4* (PBA3 (0) +PDIFF) *L4 



(III-22) 



d (HSA) /dPB3AR = -K4* (LD ( 0) +L4 *THE0TA (0) ) 



(1 11-23) 



Next the planing force partials are derived: 

HPLAN = -2*PLC3EF*RHO*ASEA*V**2*PI*THETA 
let KPLAN = (2*?LC0E?F*RH0*AREA*PI) then. 
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HPLAN 



-KPLAN *V**2*THETA 



and , 

d (HPLAN) = (d (HPLAN) /dTHET A) * dTHET A 
for this constant speed model, therefore. 



d (HPLAN) /dTHETA = -KPLAN *V (0) **2 



Plenum Gauge Pressure lift was defined to be: 



(III-24) 



HPRES = -AB*PBBAR 



(II-llc) 



as a result: 



d (HPRES) /dPBBAR = -A3 (0) 



(III-25) 



4 . Pinch Moment P arti als 



The total moments acting on the craft are: 



PTOTAL = P 3F + P3 A + PSF+PS A+ PPLAN+PPRES+PD AMP 



(II-3B) 



The 

Deginning 



partials of ir-38 will be found term 
with the forward buoyant moment: 



by 



t erm , 
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P3F = -HBF*L5 



(H-23) 



therefore: 



d(PBF) = -L5#d (HBF) 
so that: 



d(PB?) = -L5* (d (HBF) /dZ) *dZ-L5*{(d (H3F) /d THETA) * d THETA 

(III-26) 



and the partial of the after buoyant moment is: 



PB A = Lo*HBA 



(11-26) 



dP3A = Lo* dHBA 
therefore , 



d (PBA) = L6* (d (HBA) /dZ) *dZ + L6* (d (H3A) /dTHETA) *dTHETA 

(III-27) 



Seal forces partials must include the dynamic 
variation in lever arm length, which for the forward 
seal is: 



PS? = -HSF*PLSF 
or. 
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PSF 



(PBBAR *ASEAL1 ) *PLSF 



therefore , 



d (PSF) /dPBBAR = ASEAL1 (0) *PLSF (0) 



(III-29) 



and 



d (PSF) dZ = PBBAR (0) *d(PLSF*ASEAL1) /dZ 
where 



d(A3EAL1*PLSF) /dZ = PLSF (0) *d (ASEAL1) /dZ + 
ASEAL1 (0) *d (PLSF) /dZ 



ASEAL1 = WIDTH* (LD-L3*T SET A) /SIN (31) 
and 



PLSF = L3- (LD-L3 *1 HE TA) /2*SIM (31) 
therefore , 

d (PLSF ) /dZ = -1/ (2*SIN (3 1 ) 
and 

i (ASEAL 1) /dZ = WIDTH/SIff (3 1 ) 
as a result, 

d (PSF) /dZ = PBBAR (0) *( (PLSF (0) *WIDTH) /SIN (31) - 
ASEAL1 (0) /(2*SIN (31) ) 
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( r ii— 29) 



d (psf) /dTHETA = PBAR (0( *d { PLSF* AS EAL 1 ) /dTHETA 



d (PLSF* AS EAL 1) /dTHETA = 

PLSF (0) *d (ASEAL1 ) /dTHETA+ASEALI (0) *d (PLSF) /dTHETA 



d (ASEAL1) /dTHETA = - L3*ii I DTH/S IN (3 1 ) 



d (PLSF) /dTHETA = L3/ (2*SIN (31) 
as a result, 

d (PSF) /dTHETA = -P B 3 AR (0 ) * ( {PLSF (0 ) * L3 * WI DTH) /5 I N (3 1 ) - 
ASEAL1 (0) *L3/ (2*SIN (31 ) ) 

(III -30) 

For the after seal the derivatives proceed in the 
manner: 



PSA = - (PBBAR + PDIFF) *ASEAL2*PLSA 



d (PSA) /dPBBAR = -ASEAL2 (0) *PLSA (0) 



(HI -31) 

In the same manner as shown in 
paragraph, the partials of moment due to 
respect to the other variables are: 



the previous 
seals aft with 
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d (psa) /dZ = -P BB AR (0 ) * ( ( PLS A ( 0) * WIDTH) /SIN (32) + 
ASEAL2 (0) / (2*SIN (32) ) 



(III-32) 



d (PSA) /dTHETA = -PBA R (3 ) * ( (PLS A ( 0) *L4* WIDTH) /SIS ) 3 2) + 
ASEAL2 ( 0) *L4/ ( 2*S IN (32) ) 



(III-33) 



The partial of planing moment is: 



PPLAN = L7*HPLAN 



(11-30) 



d (PPLAN) = L7*d(HPLAN) 



d (PPLAN)/dTHETA = -L7*KPL AN*V (0) **2 

(III-34) 

For this constant speed model all other partlals are 
identically zero. The Plenum Gauge Pressure also has 
dynamic variation in the moment lever arm and the 
derivatives are: 

PPRE3 = -XCP*HPRES 



(IT-34) 



therefore , 



dPPRES = -HPRES (0) d (XCP) -XCP (0) *d (HPRES) 
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where 



XCP = XCPO + (LDIFF*LD)/ (2^311 BHGT) 
therefore. 



d(XCP)/dZ = LDIFF/(2*BUBHGT) 
and 



d(PPRES) = - ( (HP RES (0) * LDI FF) / (2*BUBHGT) ) *dZ 
-XCP(O) *d (HPRES) 

substituting equation (III-25) yields: 

d(PPRES) = -( (HPRES (0) *LDIFF) /(2*BUBHGI) ) *dZ 
+XCP(0) *AB (0) *dPBBAR 



Finally, the partial of pitch damping moment is: 



PDA HP = -2*PI*XLSS**2*YLSS**2*RHD*U*THETADOT/8 



let 



(H-37) 



2*PI*XLSS**2*YSS**2*RHO/8 = K D AMP 
then 



PDAMP = -KDAMP* V*THETAD3T 

as a result, since this is a constant speed model. 
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d (PDAKP) /dTHEIDT = -KDA!1P*V(0) 



( 1 11-35) 

where THETDT is the first derivative of theta with 
respect to time. 



5. Total Differentials 



Since F = M * A and the interest is in 
acceleration terms, A = F/.1 and dA = d(F)/M. the result 
being that all differentials must: be divided by craft 
mass to achieve the desired linearized heave results. 
The same general comment holds true for moment 
differentials, that is; they all must be diviied by the 
craft moment of inertia about the Y-axis. 

All differentials must be sorted by the variable 
of differentiation so that the linearized differential 
equations of motion may be written. As an example; 

dZDDOT = ( (d (HBF) /dZ) ♦ (d (HBA) /dZ) ) /MASS) *dZ + 

The State Variable representation of the system 
will be presented in a later section of this chapter, 
however it is convenient to define the State Variables 
at this time. 



t 

X = (d Z ,dZDOT , dT HE! A , dTHETDT, dMB) 



(III-36) 



and 
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Y = d P B3AR 



(III-37) 

since PBBAR is not a state variable the definition for 
dPBBAR must be substituted for dPB3AR in all equations. 



6. Total Heave Force Differentials in Z 



Here the total differentials will be found. No 
explanatory text will be included unless special 
suostitut ions requiring an aside are necessary. 

The buoyancy terms are: 



DH3Z = d(HBF)/dZ + d(HBA)/dZ 

(1 11-38) 

A term must be included which was not previously 
derived to account for the s'ensitivity of seal forces 
due to PBBAR since it is not a State Variable. 



DH5PZ = (d {HSF) /dPBBAR) * (d (PBBAR) /dZ) 
as a result. 



d (HSF) /dPBBAR = -K3* (LD (0) -L3*THETA (0) ) 



(I II- 20 ) 



and 



DHSFPZ = -K3*(LD (0) -L3*THETA(0) ) *d{PBBAR)/dZ 
similarly for the aft-er seal moment: 
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a (HSA) /dPBBAR = -K4* ( LD (3) +L4*THETA (0) ) 



( 1 11-23) 



and 



DHSAPZ = -K4* (LD (0) -t-L4*rHETA(0) ) *d (PBBAR) /dZ 
The total seal force differential in Z becomes: 

DHSZ = d(HSF)/dZ + d(HSA)/dZ +DHSFPZ + DHSAPZ 



(III-39) 

For plenum gauge pressure lift differential the 
same substitution must be made. 



dHPRES = -AB (0) *dPBBAR 
therefore. 



DHPZ = -AB (0) *d (PBBAR) /dZ 



(III-40) 



Define the sensitivity coefficient of ZDDOT to Z 

to be : 



DZZ = (DHBZ +• DHSZ + DHPZ) /MASS 



(I II -4 1 ) 



7 . Total Differential s in THETA 
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The following partials are now defined: 



DH3TH = d (HBF) /dTHETA + d (HBA) /dTHETA 



(III-42) 



DHSTH = d (HSF) /dTHETA + d (H3A) /dTHETA 



(III-43) 



DHPTH = d (HPLAN) /dTHETA 



(III -44) 

Define the sensitivity coefficient of ZDDDT due to 
THETA to be: 



DZTH = (DHBTH + DHSTH + DHPTH) /MASS 



(I 11-45) 



3- T otal Differentials due to MB 

Here since dP33AR is not a state variable the 
partial of PBBAR with respect to MB must be suostituted 
in the plenum lift and seals lift terms. 



d (HSF) /dPBBAR = -?K3* (LD (0) -L3*THETA (0) ) 



(III-20) 

oy substituting d(PB3AR|/dMB for dPBBAR The result is: 
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d (HSFP) /dMB = -K3* (LD (0* -L3*rHETA (0) ) *d (PBBAR) /dMB 



(III-46) 



similarly , 

d (HSAP) /dMB = - K4* (LD(0| + L4*THETA(0) ) *d (PB3AR) /dMB 

(1 11-47) 



d ( HP2ES) = -A3 (0) *dP3BAR 

using the previous substitution yields: 

d (HPHES) /dMB = -A3 (0) (PBBAH) /dMB 



(III-48) 

The sensitivity coefficient of heave due to MB can 
now be defined to be: 

DZM3 = (d (HSFP) /dMB +d (33 AP) /dMB + d (HPSES) /dMB) /MASS 

(III-49) 



9. Total Pitch Moment Differentials in Z 



DPBZ = d(PBF)/dZ + d(P3A)/dZ 



(III -5 0) 



The dPBBAR substitution yields the partial of 
seals with respect to pressure with respect to Z. 
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DP5FPZ = ASEAL1 (0) *PLSF (0) *d (PBBAR) /dZ 



(III-51) 



DPSAPZ = -ASEAL2 (0) *PLSA (0) *d (PBBA3) /dZ 



(III-52) 



therefore the total seal differential in Z becomes: 



DP SZ = DPSFPZ + DPSAPZ + d(PSF)/dZ + d(PSA)/dZ 



( 111 - 53 ) 



Next the Plenum Gauge pressure term is derived: 
d (PP3ES) /dZ = -HPRES (0) *LDIFF/(2*B0BHGT) 

d (PPRES) /dPSSAR = AB (0) *XC? (0) 
therefore, 

DPPRSSZ = AB (0) *XCP (0) *d (PBBAR) /dZ 
and 



DPRESZ = d (PPP.ES) /dZ + d ( P PRES Z) /dZ 

(1 11-54 ) 

The sensitivity coefficient of THETA with respect 
to Z becomes: 

DTHZ = (DP BZ + DPSZ + DPPRES)/I.YY 
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(Ill— 55) 



10 . 



DPBTH 



DPSTH 



DPPLTH 



define 



DTHTH 



to be : 



DTHDTH 



Tota l Moaent Differ en tials in THETA 



= (d (PBF) /dTHETA + d (PB A ) /dTH S TA) 



= (d (PS F) /dTHETA * d (PS A) /dTHETA 



= d (PPLAtt) /dTHETA 



the sensitivity coefficient for THETA to 



= (DPBTH + DPSTH + DPPLTH ) /III 



Also, define the sensitivity coefficient 



= (d (PD AMP) /dTHDT) /IYY 



Tot al Moaent Differentials in M3 



(I II -56) 

(III-57) 

(III-58) 

be: 

(III-59) 
in THETDT 

( 1 1 1 -6 0 ) 



56 



dPSAMB = (d (PSA) /dPBBAS) * (d (PBBAR) /dMB) 



(I II -6 1 ) 



dPSFMB = (d (PSF) /dPBBAR) * (d (PBBAR) /dMB) 



(III-62) 



DPSMB = DPSAMB + DPSFMB 



(III-63) 

Inserting the dPBBAR suds titution yields the partial of 
plenum gauge pressure moment with respect to M3: 



d (PPRES) /dPBBAR = - XCP (0 ) *DHPRBS 



d (HPRES) /dMB = -AB (0) (PBBAR) /dMB 



DPPMB = XCP(O) *AB (0) *d (PBBAR) /dMB 



( III -6 4) 

Define the sensitivity coefficient of IH2TA due to 
M3 to be: 

DTHMB = (DPSMB + DPPMB) /IYY 



(1 1 1 -6 5) 



12. Total Differentials of Air Mass 



57 



d (JIB) /dMB = (d (MB) /dPBBAR) * (d (PBBAR) /dMB) 
define the sensitivity coefficient to be: 



DMBMB = d (MB) /dMB 



(I II -6 6) 



define the sensitivity coefficient in Z to be: 



DM3Z = d(MB)/dZ 



(III-67) 



D. TOTAL MODEL SENSITIVITY COEFFICIENTS 



1 . Heave Sens it jv it y Co eff i cie nts 

The sensitivity coefficients foe the linearized 
differential equation in heave acceleration are: 

DZZ, DZTH , and DZ MB. 



2. Pitch Sensitivity. Co eff i cien ts 



The sensitivity coefficients for the Linearized 
differential equation in pitch acceleration are: 

DTHZ, DTHTH, DTHDTH, and DTHKB. 

3. P lena m ^ ir Mass Sensitivity Coefficients 
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The sensitivity coefficients for plenum air mass 

are: 



DMBZ, and DM3MB. 

4 . Ple num- Gauge Pres sure Se ns itivi ty Coef ficients 

. The sensitivity coefficients for Plenum Gauge 
Pressure are: 

DPBZ, and DPBilB. 



E. STATE SPACE REPRESENTATION 

The state matrix form of a dynamic system is: 

X DOT = AX + BU 

where X = a vector of variables representing the states 
of a sat of first order linear differential equations, 
and 0 is the vector of forcing functions, and A and B 
are the matrices of the respective gains. The form of 
the equations is: 



X1D0T = X2 



X2D0T = X3 



X3D0T = Kl*Xt + K2*X2 f K3*X3 
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for the linearized model of the CAB SES the state 
variables were defined to be : 



t 

X = (dZ, dZDOT, dTHETA, dTHETDT, di3B) 



As a result the state matrix representation of the 
linearized unforced Heave-Pitch modal is: 



X 1DOT 




0 


1 


0 


0 


0 




"xi' 


X2DOT 




DZZ 


0 


DZTH 


0 


DZI1B 




X2 


X3DOT 


= 


0 


0 


0 


1 


0 




X3 


Xh DOT 




DTHZ 


0 


DTHTH 


DTHDTH 


DTHM3 




X4 


X 5 DOT 




DHBZ 


0 


0 


0 


DM3MB 




X5 



and 



(I II -6 8) 



dPBBAB = DPBZ*X1 + D?BMB*X5 
III-69 
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IV. £I!1E fiOMAIN VALIDATION 



The linearized model developed in chapter III was 
validated against the the Six-degree-of-freedom 
nonlinear model. There were several changes made in the 
simulation techniques used, and' at the same time the 
changes in center of pressure dynamic variation and 
planing forces were entered into the 6 OOF model. 
Establishment of equilibrium for the linear model was 
accomplished off-line using a TI-59 programmable 
hand-held calculator with printer. 

The first change in simulation technique is the 
operating point at which the linear sensitivity 
coefficients were evaluated. In McIntyre's model these 
coefficients were evaluated at the initial operating 
point. In this worlc the coefficients are evaluated at 
the final steady-stare equilibrium point since this is 
cusromary in linearization work. 

The second major change in simulation was the method 
used to determine the planing force loss coefficient 
(PLCOEFF) and the establishment of the arbitrary 
location of the initial plenum center of pressure. In 
this model all heave forces are well defined (as modeled 
here) with the exception of planing force lift. As a 
result, craft weight and all heave forces except that 
due to planing lift are summed and the residual assumed 
to be due to the planing phemonema. This residual is 
then used to determine (off-line) the value of PLCOEFF. 
Note that this method also accounts for any lift due to 
aerodynamic forces on the bow that were disregarded. 
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Next all moments are summed about the C.G. except the 
moment due to plenum pressure lift. The residual moment 
is assumed to be due to plenum pressure. This residual 
is used (again off-line) to establish the initial center 
of pressure. This accounts for any moments not 
previously accounted for (such as aerodynamic forces on 
the bow.) McIntyre used the planing lift coefficient to 
force the model into agreement with the 6 DUF model for 
natural pitch frequency and used a large, arbitrary, and 
static moment to force the models into steady state 
pitch angle agreement. It was felt that using the 
present method of determination would leave the 
equations in a form that would be more amenable to 
extraction of the factors most directly affecting craft 
characteristics as modeled herein. 



A. LINEAR PROGRAM DESCRIPTION 



User's instructions have been included as appendix A. 
The program flow can be isolated into four main 

sections. 

1. First the values for craft weight, draft, pitch 
angle, plenum pressure and speed are entered. To 
initialize plenum air mass, draft is used to solve for 
initial plenum air volume, which is then used with 
initial plenum gauge pressure to solve equation 11-15 
for plenum air mass. Equation 11-12 is then solved for 
Qin. Qout is set equal to Qin and equation 11-13 is 
solved for seal leakage area. All other heave and pitch 
nonlinear equations are solved and their values are 
listed along with residuals of heave aid pitch 
summations. These residuals are added to the initial 
values of C.G. and angular acceleration and tne initial 
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values are saved for the time domain simulation section. 

2. The same operations are carried out for tae steady 
state values of weight, draft, pitch angle, and plenum 
pressure. The residuals, forces and moments are again 
listed. 

3. These steady state values are used to evaluate 
the sensitivity coefficients developed in chapter III 
which are also listed. 

4. Finally the state variable solution is simulated 
using a four-point fixed step Runge-Xutta algorithm and 
listed for the print interval desired. The option is 
included to output either printer plot grapnicaL display 
using the NPS printer plot subroutine PLOTP or the X-Y 
Versatec plot subroutine DR ASP . The latter should be 
used for final runs only, with PLOT? used on 
intermediate runs. Provisions are made for up to 300 
discrete output points for the time domain tabular 
listing . 



3. SIX-D EGREE-OF-FREEDOM MODEL SOLUTION 



Two changes were made in the 6 DOF model to account in 
a more analytic manner for the observed rendancy of the 
XR-3 to pitch down as speed increased. 

The first change was to remove an empirical equation 
used to correct the center of pressure L5] and replace 
it with the dynamic center of pressure variation 
developed in chapter II. 

The second change * as to include in the sidewall 
subroutine the nonlinear equation for planing forces, 
wicn the loss coefficient adjusted to make up (as 
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closely as possible) foe lift loss due to the revised 
modeling of the plenum chamber. 

Finally/ the 6 DOF model was ran at both 6722 lbs and 
6050 lbs weight to establish the initial and final 
operating points, followed by the weight removal 
transient run for tine domain validation of the linear 
heave-pitch modal. 

It should be noted that further study of the 6 DOF 
model is needed to evaluate the results of the two 
changes and to more closely establish the arbitrary 
initial location of the center of pressure. The large 
amount of CPO and turn-around program times required to 
run the nonlinear model precluded including that study 
in this work. Also, for the nonlinear model, 
X-direction thrust was held constant during the run with 
the result that speed was 30.33 knots and still 
increasing at the end of the five second run. This 
effect will be further discussed in the next section. 



C. COMPARISON OF H2SULT5 

The graphical output for the five second runs are 
shown in figures 5 through 8 for the linear model and 
figures 9 through 12 for the nonlinear model. The fixed 
integration step size was 0.005 sec for the linear model 
and 0.001 sec for the nonlinear modal. 

1 . Plan urn Gau<je Pressure Tr ansi ent 
(see fig. 5 and 9) 
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In general both models exhibited a rapid initial 
redaction in pressure followed by a smooth transient 
back towards the final operating pressure. For the 
linear model the minimum occured at 100 milli-seconds, 
and 110 milli-seconds for the 5 DOF model, with the 
respective pressures being 21.38 and 21.30. The final 
values were 24.31 and 24.08 for the linear and nonlinear 
model respectively. 

As can be seen from the graphical data excellent 
agreement was acheived between the two models. 

2 . Dra ft Tr ansi en t 

(see fig. 6 and 10) 

Both models exhibit a smooth exponential type 
transient towards the final, lower steady-state draft. 
The linear model transient settles towards a lower value 
(6.04 inches as compared to 6.44) which is in part 
attributed to a difference in modeling of the seal 
forces. The nonlinear model uses a complex combination 
of hydrodynamic and hydrostatic lift where the linear 
model uses a simple pressure differential effect. Also, 
since the nonlinear model is settling to a significantly 
lower pitch angle (see para. 4 below) , the planing force 
effect is reduced in that model. 



3. C.G Acceleration 



(see fig. 7 and 11) 

Both models show a rapid positive transient in 
acceleration that closely follows the negative transient 
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in plenum pressure. Both models reached the initial 
peak positive C.G. acceleration at approximately 20 
milli-seconds after the plenum pressure completed the 
negative transient. The linear model appears to have 
greater coupling between pitch and heave looking at the 
C.G. acceleration, but this is not actually supportable 
since the linear model has greater pitch amplitude 
excursions about the exponential pitch transient to the 
lower steady state value (see fig. 8 and 12.) In 
general the two models show excellent agreement in 
transient and damped C. 3. acceleration behavior. 

4 . Pit ch Transient 

(see fig. 8 and 12) 

It is in this area that the two models are in 
greatest disagreement. Both exhibit a general damped 
sinusoidal behavior about a center exponential transient 
to a lower steady state value. In both cases the 
initial excursion is downward.' The natural pitch 
fequency of linear model is approximately 4.3 radians 
per second, whereas the nonlinear model exhibits a 
natural pitch frequency of 6.98 radians per second. 
There is felt to be one major factor contributing to 
this large discrepancy, which is the inclusion of the 
planing force in the nonlinear model. Data presented by 
Thaler and Gerba 0 4} showed a natural pitch freguency of 
5.8 radians per second and much higher oscillations 
about the exponential transient. It is felt that the 
simplified planing forces of the linear model are 
included in the much more complex hydrodynamic aquations 
of the nonlinear model and are therefore redundant in 
that modal. Due to time constraints, that avenue of 
investigation is left for future studies. The linear 
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model is tending towards a steady state value on the 
order of 0.45 degrees whereas the nonlinear model is 
tending towards approximately 0.30 degrees. This large 
discrepancy is partially due to the redundant planing 
force in the nonlinear model and the fact that speed is 
also increasing in the nonlinear model which tends to 
accentuate the planing force and force the craft to 
pitch down. 

In the next chapter a signal flow graph will be 
developed and Mason's Gain Rule applied to explain the 
linear model characteristic response in a more 
analytical fashion. 
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V. 



SIGNAL FLOW’ GRAPH AMD MASON'S GAIN ROLE 



A. INTRODUCTION 



The State Variables defined in Chapt II are used to 
develop a signal flow graph of the craft heave-pitch 
dynamics. Mason's Gain rule is applied to the signal 
flow graph to derive analytical S-domain polynomials 
expressing the transfer functions for pitch moment into 
pitch angle, heave acceleration into draft and the two 
cross gains. Using typical gains calculated in the time 
domain program developed in Chapt III the roots of these 
polynomials are located and it is shown that pitch 
dynamics may be vary closely approximated by a 
second-order system. Heave dynamics are also 
simplified. 

Using a simplified approach to the frequency response 
of the craft to non calm-water conditions, the basics 
are developed for later sea-state model work. A dynamic 
system modeling program available at tae Naval 
Postgraduate School is used to generate Bode plots for 
the craft's frequency response in Heave and Pitch under 
the simplifying assumptions. 

B. SIGNAL FLOW GRAPH ANALYSIS 

Fig 13 was generated using the sensitivity 
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coefficients defined in Chapter III. From that graph the 
following feed-forward and feedback loops and paths were 
identified: 

1. Feed Forward, U(THETA) into THETA. 



PI = (1/IYY) * (VS**2) 

2. Feed Forward, Q(Z) into Z. 

P2 = (1/HASS) * ( 1/S**2) 

3. Feed Forward, L T (THETA) into Z. 

P 3 = (1/IYY) *DZT H* ( 1 /S**4) 

4. Feed Forward, U(Z) into THETA. 

P4 = ( 1/HASS) *DIHZ* ( 1/S**4) 

5. Feedback Loops. 

LI = DTETH/S**2 



L2 = DT HD TH/S 



L3 = (DZIH) * (DTHZ) /S**4 



L4 = (DZMB) * (DMBZ) /S **3 



\ 
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L5 = DZZ/S**2 



L6 = DMBMB/S 



L 7 = (DZTH)*(DMBZ) * (DT3MB) /S**5 



C- MASON * S SAIN ROLE APPLIED 



Mason's gain rule for a closed multi-loop system is: 



I = SUM (Pk*DELTAk) /DELTA 
where , 

T = transmission gain from an input to an output 
node. 

DELTA = the graph determinant 
DELTAlc = cofactor of the kth path 

Pk = gain of the direct path from input to output 

Using Mason's gain rale the graph determinant (DELTA) 
is found to be of the form: 

D (s) = (35*S**5 + B4*s**4 + 33*s**3 + B2*s**2 + B1*s ♦ 

BO) /S**5 

The forward direct path U (THETA) to THETA is: 

N 1 (s) = (1/IYY) * (1/S**5) * (A3*s**3 + A2 *s**2 + A1*S + AO) 

The forward path U (Z) to Z was found to be: 
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N2(s) *= (1/MASS) *(1/S**5) *(A3*S**3 *■ A2*s**2 + A1*s + 

AO) 

The forward direct path O(THETA) to Z is: 

N3(S) = ( 1/IYY) *DZIH*(S + A0)/S#*4 

The forward path for draft coupling into pitch { U ( Z) 
to THETA is: 

N 4 = (1/MASS) *DTHZ* (S *■ A0)/S**4 



D. SOLUTION OF THE INDIVIDUAL POLYNOMIALS 



The gains used below were taken from a typical time 
domain run as developed in chapter IV. 



DZZ = -2178.0 
DZTH = -121.3 
DZMB = -5499. 0 
DTHZ = 28.78 
DTHTH = -22. 3 1 
DTHMB = 72.1 
DMBZ = -28.79 
DMBMB = -73.11 



1. Solution of the Graph Determinant 



Inserting the numerical values for the typical 
sensitivity coefficients listed above and using only the 
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dominant terms where possible, the roots of the graph 
determinant were found to be: 

DELTA = (S**2 + 71.8*S * 2066)*(S + 0. 5 086) * (S* * 2 + 
0.72*S +22.74) /S**5 

These factors represent a pair of compLex poles 
located at -35.71+ j27.34 , a real pole located at 
-0.5085 , a pair of complex conjugate roots located at 
-0.360 + j 4.755 , and a fifth order zero located at the 
origin . 



2. S olution of Ml (S ) 

N 1 (S) = (1/ITY)*(S + 0. 427)*(S**2 + 72.68*S + 2146) /S**5 

This polynomial factors into a real zero located 
at -0.427 , a pair of complex conjugate zeros located at 
-36.34 + j 28.74 , and a fifth order pole located at the 
origin. 



3. S olution of M2iSl 

N 2 (S) = (1/MASS) * (3**2 + 0.7721*S + 22.35)*(S + 

72.38) /S** 5 

Which factors into a pair of complex conjugate 
zeros located at -0.3856 + j 4.712 , a real zero at 

-72.38 , and a fifth order pole located at the origin. 
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4. Solutio n of N3 (S ) 



N3(S) = (1/IYY) *DZTH*(S ♦ 73 . 1 1 ) /S**4 



Which represents a. real zero at -73.11, and a 
fourth order pole at the origin. 



5. Solu t ion of N4 (S ) 



N4(S) = ( 1/MASS) *DTHZ* (S > 73 . 1 1 ) /S**4 

Which has the same pole and zero locations as 

N3 (S) . 

E. SIMPLIFICATION OF TRANSFER FUNCTIONS 



1 • Tit SL 



T1 (S) = THETA (S) /U (THETA, S') 

Therefore from Mason's Gain Rule: 

T 1 (S) = N1 (S) /DELTA 

T1 (S) = (1/IYY) * (S + 35 . 34 + 428. 7) * (S + 0.427)/ 

(S + 35.1 + j 27 . 84) * (S + 0 .508) * (S * 0.36 ± j 4.755) 
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To simplify let tie complex conjugate zeros at 
-36.34 + j28.7 cancel the poles at -35.9 + j27.84 , and 
the real zero at -0.427 cancell the pole at -0.508. 
The result is the simplified second order system: 

T1 (S) = (1/IYY) / (S**2 + 0.72*s + 22.74) 



Inspection of the time domain pitch response with 
respect to natural frequeny and damping confirms this 
basic underdamped second order response. 



2. T2 (S ) 



T2(S) = Z (S)/0 (Z,3) 



12 (S) = (1/MASS) * (S + 0.3856 + j4.712)*(S + 72. 38)/ 

(S** 2 +35.1 + j 2 7.34) * (s + 0.5085)*((S +0.3S + 
j4 . 755) 



Canceling the pole at -0.36 + j 4.755 with the 
zero at -0.3856 + j4.712 results in the simplified 

response: 

r2(S) = ( 1 /M AS S) * (5 + 72 . 3 3) / (S + 0.508)*(5**2 ► 71.8*S 

+ 2066) 



Further simplification can be 
72.33) /(S + 0.508) is approximated by 

(which is valid for longer time period 
As a result T2(S) becomes: 



made if (S + 
72.33/(3 + 0.508) 
considerations. ) 
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T2 (S) = (72. 33/MASS) /(5*0. 508) * (S**2 ♦ 71.8*S + 2066) 

which has poles located at -0.508, and -35. 9* j27.83 

This is in quite close agreement with the results 
obtained by Thaler and Serba 14 for the same airflow 
condition (Kq = 1 .0) for rheir heave-only Linearized 

model. Their re-sults were a real pole at -0.33 and a 
complex pair at -35.1* j34.4 . This is of the form: 

1(5) = K/ (S + 0 . 38) * (5**2 + 70.2*S + 2415) 

3. T3JS1 



T3 (S) = Z (S) /U (THETA) 



Here it is noted that there are no pole zero pairs 
that conveniently cancel so that the transfer function 
becomes : 



T3(S) = (1/IYY) *DZTH*S* (S + 72.33) / 

(S + 35.1 + j 27.84)*(s + 0 . 5085) * (S + 0.36 * j4. 775) 

Making the simplification used in rhe previous 
section, the result is: 

T3 (5) = (72. 33/IYY) *DZTH*S/ (3**2 * 71.8*S + 2055)* 

(S + 0. 508) * (S** 2 + 0 . 72* S + 22.74) 

Which can be further simplified by looking at the 
long term response vice initial few milli-seconds to 
(S/(S * 0.508)) = 2 so that: 
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T3(S) = (142/IYY) *DZTH/ (s**2 + 71.8*3 + 2056)*(3**2 + 
0 . 7 2*s ♦ 22.74) 



4. Tt (S) 

Since the transfer fundtions are the same except 
for the gains the solutions are the same except for the 
gains : 

T4<S) = ( 1/MASS) *DTHZ * S* (S *■ 72. 33) /DELTA 

where the same simplifying assumptions may be applied. 



?. SIMPLIFIED BODE PLOT ANALYSIS 



t. Sim p lify ing Assu mption 
(see fig. 12) 

The assumption is made that the craft linearized 
frequency response transfer functions can be separated 
into two transfer functions. 

The first block G1(W,U,...) expresses the 
disturbing moments and accelerations generated by the 
craft interacting with the waves. This function was 
developed by Thaler and Gerba C14J for the heave-only 
model- In that simpler model the function G1 was a 
complex function of encounter frequency, angle of 
incidence, craft length compared to encounter wavelength 
etc.. This portion of the overall frequency response of 
phis model is left for future work due to time 
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constraints. 



The second block 32 (We) represents the heave-pitch 
dynamics of the model as developed in chapters II, III, 
and IV, which respond to the disturbing moments and 
accelerations from block G1. 

The Transfer functions T1 and T2 (unsimplif ied) 
were simulated using the IBM-360 simulation language DSL 
with S = jW with the following results. 

2. Bod e Plots for Pitch Dynamics 
(see figs. 15 and 16 ) 

Figure 15 represents the normalized (numerator 
gain constants set to 1) gain in dB's versus Omega. 
Mote that the abscissa scale is Omega in powers of 10 
since this is a convenient means to achieve a 
logarithmic scale on the X-Y plotter output. is can be 
seen the response strongly supports the simplification 
to second-order system. These plots represent the full 
transfer function T1 (jW) vice the simplified system. 
Figure 16 further supports the simplification as the 
phase plot starts at zero degrees and shows only a very 
slight positive increase prior to the rapid change 
toward -130 characteristic of an underdamped system. 
Also note that above the damped natural frequency the 
magnitude plot falls at -12 dB per octave. 

3. Bod e Plots for Heave Dynamics 
(see figs . 17 and 1 3) 
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Figure 17 represents the normalized magnitude 
versus Omega for the heave dynamics. The system is 
characteristic of an overdamped third order system with 
small gain increases at the craft's natural pitch 
frequency, which represents the pitch dynamics feeding 
back into heave dynamics. However the gain magnitude 
curve has already fallen 20 dB below it's lower 

frequency response gain before the pitch coupling begins 
to manifest itself. As can be seen in fig. 18 the phase 
plot shows the same third order pole response with 
slight second order zero influence at the natural pitch 
f requency . 
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VI 



. CONCLUSIONS AND RECOMMENDATIONS 



As shown in Chapters IV and V the simplified 
Pitch-Heave model developed herein is a very reasonable 
simulation of the XR-3 test craft for the given 
operating conditions. Heave, Pressure, and C.G. 
Acceleration were in very close agreement. Pitoh on the 
other hand showed some difference in natural frequency 
and damping factor. It is felt that this difference is 
primarily attributable to the inclusion of a redundant 
planing force calculation in the 6 DOF model. 

Historically the results of the 6 DOF computations 14 
for pitch response have been very similar to the results 
for the linear model, in both natural frequency and 
damping characteristics. This can be attributed to the 
fact that the additional planing force did not exist in 
the original 6 DOF program. 



Since one of the major concerns with the 6 DOF model 
is the excessive computational times involved in 
sea-state simulations the conclusion may be drawn that 
eventual substitution of the sidewall and seal force 
linear model equation in the 5 DOF program will 
substantially reduce CPU time, with some resultant loss 
in computational accuracy. By way of comparison, 
disregarding compiling time and output time, the linear 
model required approximately 9 seconds of CPU time for 
the 5 second simulation while the 6 DOF model required 9 
minutes of CPU for the 5 second run. This is a 
computation time improvement on the order of 60 times. 

From the time-domain plenum gauge pressure data it can 
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be seen that the dynamic variation in effective plenum 
roof area has a noticeable effect on the plenum pressure 
transient behavior as compared to previous results 14 . 
The initial pressure down transient was faster, dropped 
to a lower value , and began it’s recovery more quickly. 
Therefore the plenum roof dynamic variations strongly 
affects the C.G. acceleration characte ristics of the 
craft and therefore it's transient behavior. 

The dynamic center of pressure variation developed 
herein proved to be a very draft sensitive phenomenon. 
Although the length o£ the other lever arms as modeled 
herein change by a similar amount, percentage of 
variation is radically different. For the seals the 
variation is on the order of inches in 10 feet, where 
the variation in lever arm for plenum pressure lift is 
on the order of inches in a single foot (or lass.) This 
implies that dynamic G.P, variation is the dominant 
factor in steady state pitch angle versus steady state 
draft. 

Planing forces as included in this model are now in a 
form where the basic factors as far as craft 
construction characteristics are readily recognized. 
This force provides part of the answer to craft 
pitch-down characteristics as initially investigated by 
Reidel C93 , with the remainder of the moment being 
developed by the C.P. variation for pressure. The 
planing force provides a clue as to how, during 
construction, the craft’s pitch amplitude excursion 
characteristics may be controlled. 

Several areas for future investigation have become 
apparent in this work. 

1. Verification of the effect of dynamic center of 
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pressure variation and planing forces as incorporated in 
the 6 DOF model over the entire speed operating range. 

2. Incorporate the simple modeling for sidewalls and 
seals into the 6 DOF model and investigate the trade-off 
between CPU time and computational accuracy. 

3. A closer look at the seal forces used herein. It 
is possible that a simple gain factor needs to be 
included in these equations in order to optimize the 
steady-state seal forces with respect to actual load 
data as measured during experimental test runs of the 
craft 10 . 

4. Incorporate sea state forces and moments into the 
simplified model by simulation of the forces and moments 
that are generated by craft interaction with the sea. 
Evaluate the accuracy of the method used by comparison 
with data generated by the simple model and the 6 DOF 
model in sea-state conditions. 

5. Analytical investigation of the craft physical 
characteristics that have predominant influence on the 
natural pitch frequency and damping by using the gain 
terms from the results of the signal flow graph analysis 
done in Chapter V. 

6. Use the simplified model to develop a ride-comfort 
control system. lest the resulting design on the 6 DOF 
model, and compare the results obtained with those 
suggested by 30 HR C73 . It has been known empirically 
for some time that the craft sxhioited simple second 
order behavior in the pitch mode. Mow that this fact 
has been analytically verified the implications for an 
automatic control system for pitch response are 
powerful. A simple oontrol system can be visualized 
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that does not involve either complex optimal control 
system design work or hardware. A simple acceleration 
and rate feedback control system is implied, with the 
possibility of adding a speed input to the system to 
compensate for the localized linearity of the craft 
characteristics. 



32 



APPENDIX A 



LINEAR 5 E S MODEL USER'S GUIDE 



A. General Comments 



The program as listed following Appendix 3 is intended 
for useage via the hot card reader at NPS. 

It is intended to run in Fortran G with Plotting 
option (PROC = FOP.TCLGP . ) 

The program is structured on the assum-ption that the 
initial and final state variables are in equilibrium 
prior to data entry, but will run regardless. The force 
and moment residual listings will indicate the value of 
imbalance, if any. 

The craft physical dimensions are alterable by 
changing the cards in the craft dimension section of the 
main program, which are keyed to the craft geometry 
herein. 

3. Data Deck Setup. 

The /* and //GO.SYSIN DD * cards are inserted 
following the and card in Subroutine VERSAP, foLlowed by 
the data deck as described 

I.Plot option card (IPLOT) , format II. iPLOr = 1 = 

VERSATEC X-Y plotter output. 

IPLOT = 0 = PET PLOT output 
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2. List option card (ILIST), Format II 

ILIST = 1 = Forces and moments, with residuals and all 
sensitivity coefficients listed. 

ILIST = 0 = Input craft operating conditions only 
listed . 

3. Time and Integration variables. 

TI, TF, NPTS, NSTEP Format (2F1 0. 4, 2110) 

Ti = problem start time (normally 0.0) 

TF = Problem stop time (seconds) 

NPTS = number of tabular output points desired. 

NSTEP = number of integration steps. 

note: it is best to make NSTEP an integer multiple 

of NPTS, and never less than NPTS. 

4. Planing coefficient and arbitrary center of pressure 
input, i. e. HEAD ( ) AKP, PLCOEFF 

Format = 2F10.4 

5. Initial Conditions card. 



The initial craft conditions are 
following order: w eight (Lbs .) , Draft 

(Deg.)), PEBAR (psf) , and speed (fps.) 



read in in the 
(in) , Pitch angle 
Format is 5F10.4 



6. Final Conditions. 



The final conditions are read in the following order 
and format: weight (Lbs), Draft (in). Pitch angle 

(Deg), and PBBAR (psf) , Format 4F10.4. Since this is a 
constant speed model, initial and final speed are 
assumed equal. 



7. Label Cards. 



Label cards for VERSAIEC output are the last 10 cards. 
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There are five graphs in this mode. Draft(in) , C. G. 
Acceleration (G.' s) , Pitch Angle (Deg.), Pitch 
Acceleration (Deg/SEC**2) , and PBBAH (Each graph 
requires two label cards. Format (6A8) , which print on 
two lines below the graph. The graphs are output in the 
order listed above. 






85 



APPENDIX B 



DRAWINGS AND FIGURES 
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Figure 2 - ASSUMED GEOMETRY <C0MT> 
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Figure 3 - ASSUMED GEOMETRY <CONT> 
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Figure 4 - PLENUM CHAMBER MODELING 
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figure 5 - LINEAR PLENUM PRESSURE TRANSIENT 
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Figure 6 - LINEAR DRAFT TRANSIENT 
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Figure 7 - LINEAR C. G. ACCELERATION 
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Figure 8 - LINEAR PITCH TRANSIENT 
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Figure 10 - 6 DOF Dfi APT TRANSIENT 
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Figure 11 - 6 DOF C.G. ACCELERATION 
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Figure 12 - 6 DOF PITCH TRANSIENT 
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